Ferromagnetic and antiferromagnetic spin fluctuations and superconductivity in the 

hep-phase of Fe. 
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High purity iron, which transforms into the hep phase 
under pressure, has recently been reported to be. super- 
conducing in the pressure range 150-300 kBarfcl. The 
electronic structure and the electron-phonon coupling 
(X p h) are calculated for hep iron at different volumes. A 
parameter-free theory for calculating the coupling con- 
stants A s / from ferromagnetic (FM) and antiferromag- 
netic (AFM) spin fluctuations is developed. The calcu- 
lated A s / are sufficiently large to explain superconductiv- 
ity especially from FM fluctuations. The results indicate 
that superconductivity mediated by spin fluctuations is 
more likely than from electron-phonon interaction. 



Recent observations of superconductivity-, in several 
magnetic or nearly magnetic materialsou have re- 
newed the interest, in theories of paramagnon mediated 
superconductivitycTD. Progress in metallurgy has been 
very important for these discoveries, since this type of 
superconductivity is sensitive to perturbations and im- 
purities. In particular, the recent report of superconduc- 
tivitvJn the hep phase of iron, at a pressure of about 150 
kBaitH is interesting, since it is a non-compound material 
where high purity might be easier to achieve. A sub- 
sequent theoretical work concluded that this supercon- 
ductivity could be due to conventional electron-phonon 
coupling, but that the disappearance of superconductiv- 
ity at .higher pressure can not be explained by the same 
meanstl. Spin fluctuations, ferromagnetic (FM) or anti- 
ferromagnetic (AFM) have an important role to suppress 
superconductivity in this scenario. The alternative expla- 
nation, where superconducting pairing is caused by spin 
fluctuations, has been put forwardJnr nearly magnetic 
materials such as ZrZn.2 or SrRuOaErll Here we investi- 
gate the possibilities for superconductivity mediated by 
spin fluctuations in hep iron. 

The electronic structure of bec and hep Fe has been cal- 
culated using the linear muffin-tin orbital-method in the 
local spin density approximation (LSDA)oEI Non-local 
potential corrections included in the generalized gradient 
correction (GGAHl) scheme is known to be crucial for a 
correct description of the fee vs. bec stability of Fe, but 
other properties are not very [different in the two types of 
density functional potentials^. Since we are not focus- 
ing on structural stability in this work we use the LSDA 



potential. As in refH it is found that two AFM configu- 
rations are more stable than nonmagnetic or FM config- 
urations of hep Fe at high pressure. The first (AFM-I) 
has opposite polarization on each alternate z-plane layers, 
while the second (AFM-II) has opposite polarization on 
each layer perpendicular to the x-axis. The AFM-I con- 
figuration can be described within the normal hep unit 
cell, where the bands are calculated in 252 irreducible k- 
points. The AFM-II configuration is described in an or- 
thorhombic unit cell containing 4 atoms, and the bands 
are calculated in 200 k-points within 1/8 of the Brillouin 
zone. The basis set includes states up to £=3. The c/a- 
ratio is fixed to 1.59 for all volumes. This is close to the 
calculated average value within a wide pressure rangefij. 
The d-band is the dominant character at the Fermi en- 
ergy (Ep) and it is involved in the magnetic orderings. In 
the case of spin fluctuations we apply staggered magnetic 
fields locally on each site. The calculations for the three 
types of magnetic configurations are made independently 
of each other. 

The AFM-I configuration develops stable moments for 
the lattice constant a > 4.81 a.u., and AFM-II for a > 
4.61 a.u., approximately. The total energy of the AFM-II 
state is nearly 2 mRy/atom lower then the AFM-I con- 
figuration at the largest, volume. These results are in 
fair agreement with ref.Ej despite the use of different po- 
tentials. The calculated minimum in total energy of the 
bec structure is found at lattice constant which is ~ 3 
percent smaller than the experimental one. The bulk 
modulus calculated at the experimental lattice constant, 
about 1.5 Mbar, agree well with experiment, while it is 
about 2.2 Mbar at the calculated equilibrium. Such er- 
rors are common when using LSDA for 3d metals. The 
present results give a lower total energy of the hep phase 
(compared to bec) for a < 4.75. 

The electron-phonon coupling parameter X ep is calcu- 
lated from the electronic structure; 



A ep = NV 2 /Mlj 2 = NV 2 /K 



(1) 



where N is the density of states at Ep, M the atomic 
mass, to a phonon frequency and K a force constant. It 
is related to the change of total energy E due to an atomic 
displacement r. K = d 2 E/dr 2 . The Hopfield parameter 
NV 2 contains V, the Fermi surface average of the matrix 
element < ip k (dv(r) / dr)ipk >, the first order change in 
potential due to the displacement, where ipk is the wave- 
function at k. For elementary metals it can to a good 
approximation be calculated within the rigid muffin-tin 
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approximatiorJll. This matrix element could also be ob- 
tained from the change in band energies, when the 
structure is deformed as in a 'frozen' phonon. The force 
constant is fitted to the experimental value of the Debye 
temperature 9d for bcc Fe, and then scaled ta_hcp and 
other volumes by the calculated bulk modulusE£l. 

For FM or AFM spin fluctuations we make a develop- 
ment analogous to the case of frozen phonon calculations. 
Instead of applying a force on an atom to obtain a dis- 
placement as in frozen phonon calculations, we apply a 
magnetic field £ to obtain a magnetic moment m. This 
analogy can be extended to the calculation of mass en- 
hancement due to spin fluctuations, A s /. 

A s/ = NV 2 /K — N < ijj k {dv{r)/dm)ip k > 2 /{d 2 E/dm 2 ) 

(2) 

In practice, the calculations are made at two different 
configurations, one with an applied field and one with- 
out, to give two band structures and two magnetisations. 
If the changes of the free energy (E) are harmonic, i.e. 
quadratic functions of the induced magnetic moment, 
we have E^ = Eq + Jm 2 where to is the moment (in- 
duced by the field £) per atom. E^ is the free energy 
for the configuration with moment to, which is the total 
energy plus the term moment times field. The constant 
J = \d 2 E/dm 2 — (E^ — E )/m 2 is calculated from spin 
polarized results, where the applied fields range from 1 to 
10 mRy. The local Stoner enhancement, S, is defined as 
(ej — eo)/£, where is the increase of the exchange split- 
ting at an atom (obtained from the logarithmic derivative 
of the d-band) induced by the field £. The ^-values are 
generally smaller for large £, indicating some aversion 
against large moments. Non-harmonic variations of J as 
function of field are found, especially for FM cases. This 
tendency is less pronounced for AFM configurations. 

The matrix element is estimated from the change in 
band energies < ^kdv^pu > 2 = (e™— e°) 2 , where is the 
band energy at point k for the configuration belonging to 
to. In the harmonic limit we obtain; 



A s/ = N < (e™ - e° k ) 2 > FS /2(E m - E ) 



(3) 



Here < >fs means a Fermi surface (FS) average. 

The same formalism can be used for FM and AFM 
fluctuations. An alternative determination of A s / for the 
FM case from the Stoner model and paramagnetic band 
resultaiS gives A s / = ^SS 2 . The usual Stoner parame- 
ters are related by S = (1 — S)' 1 , and S — NI, where / is 
the exchange integral. The two methods give very similar 
values for X s f as long as the Stoner enhancement is the 
same in paramagnetic and spin polarized calculations. 
However, it is noted that the spin polarized results give 
somewhat larger Stoner enhancements than the param- 
agnetic calculations. _ 

Our calculated X ep are of the same order as in reffi, but 
the decrease as function of pressure (P) is faster. This 



is mainly caused by our P-dependence of the force con- 
stant K. The exchange integral I, used for the FM Stoner 
factor, calculated from paramagnetic results, varies from 
0.065 to 0.067 Ry/atom from^he largest to the smallest 
lattice constant. Mazin et. aB obtained 0.075 Ry/atom 
in their fixed spin method using GGA. This together 



with a slightly smaller 
tors smaller than in ref. 



OS makes our FM Stoner fac- 
, but it makes no essential dif- 
ference for the conclusion regarding superconductivity 
based on electron-phonon coupling. In particular, our 
calculation of A s / without parameters penpaits to de- 
termine the parameter a introduced in ref.B to fit the 
transition temperature near the volume corresponding to 
a=4.76 a.u.. With A s/ = aS, a=0.13 and S=4.6, they 
obtain X s f — 0.6, while we have A s / = 0.54 at this vol- 
ume, cf. table 1. The overall fair agreement between the 
two methods of estimating X ep and A s / (the latter tends 
to suppress superconductivity) leads to similar conclu- 
sions for electron-phonon based superconductivity as in 
refJj. Calculated electron-phonon coupling in nonmag- 
netic simple cubic Fe is of the same order and supercon- 
ductivity of the order of 10 K has been predicted!!!. 

We now turn to the possibility that superconductivity 
is based on spin fluctuations. A comparison of the lo- 
cal Stoner enhancements for FM and AFM cases, reveals 
that the enhancements can be much larger in the latter 
cases, in particular for the volumes near AFM instabil- 
ities (cf. Table 1-3). The AFM enhancements decrease 
rapidly for larger or smaller volumes. The enhancements 
for hep Fe are calculated within a wide range of volumes, 
although it should be noted that bcc is stable structure 
at the largest volumes. The decreasing Stoner enhance- 
ments at large volumes, within the AFM stability region, 
indicate that it is relatively difficult to increase the lo- 
cal moments further by applying fields. The same effect 
is probably behind the nonlinear behavior at small vol- 
umes, and for FM cases at all volumes, where there is a 
saturation of the moments for the largest applied fields. 
The only exception to this behavior at larger fields, is 
for AFM fluctuations for volumes just smaller than the 
critical volume for the AFM instability, where the non- 
linearity is reversed. The Stoner enhancements are then 
larger for the largest field, and J are 'softer'. It is as if 
a soft AFM mode will approach 'zero' field regime near 
the critical volume for the AFM instability. The rela- 
tively small FM enhancements reflect the fact that hep 
Fe is never near a FM instability in this range of volumes, 
but a metastable FM state is found at large volumeO. 

Despite much larger S for some AFM cases, there are 
no large differences between the A s / in the FM and the 
two AFM cases. One explanation is that the matrix el- 
ement for FM fluctuations is more efficient, because as 
is easy to imagine, a FM field applied in a system where 
one band is dominant at Ep, will split all bands about 
equally. The same difference (e™ — £°) will appear al- 
most everywhere over the FS. For induced AFM config- 
urations it is more difficult to visualize these differences. 
Positive and negative differences arc likely, which means 
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that small differences should exist at some sections of the 
FS and make the averaged matrix element smaller. 

The energy of the fluctuations, Huj s f, is estimated 
from a Heisenberg model, where the exchange integral 
corresponds to our parameter J. The k-dispersion is 
lu = w s /(l — cos(ka)) for FM, and cj = uj s fsin{ka) for 
AFM fluctuationsEJ. The characteristic u s f to use in an 
estimate of the superconducting transition temperature 
should be some k-point average. For ferromagnetic fluc- 
tuations it has been proposed that fko B f = l/(4NS)u. 
This corresponds to J/4, since m 2 /(NS), the estimated 
free energy from the Stoner model, is equal to Jm 2 in 
the present approach. 

For calculating the superconducting transition temper- 
ature we use a weak, coupling BCS-like formula for para- 
magnon couplingLJu. 

k B T c = (hu.f/1.2)exp(-Q. + X ep + A 5 /)/A s/ ) (4) 

The results using this formula are only approximate, but 
it turns out that they are of reasonable amplitude. The 
results for the three types of fluctuations are shown in 
the figure. The Tc as function of lattice constant in the 
case of FM fluctuations has a triangular shape with a 
maximum reaching 25 K at the largest volume. How- 
ever, the stable structure is bcc for a larger than about 
4.75 a.u. according to the calculations. Experimentally, 
the transition pressure is about 150 kBar, which would 
put the transition further away from the bcc equilibrium, 
somewhere towards 4.6 a.u.. From the calculated results 
one expects a Tc of the order 12 K at the hep transition, 
while it levels off rather slowly at higher pressure. The 
difference in P between a=4.7 and 4.4 a.u. is of the order 
700 kBar, while Tc at the smallest lattice constant still 
is about 1 K. According to experiment it is only within 
a pressure range of 150 kBar that superconductivity is 
observed. 

Superconductivity from AFM fluctuations appears 
only for volumes (on the nonmagnetic side) near the 
AFM instability. The Tc's are lower and the pressure 
range narrower than in the FM case. The AFM-II con- 
figuration extends into smaller volumes than AFM-I, fur- 
ther away from the bec-hep transition. The drop in 
Tc as function of pressure is more significant and corre- 
sponds better to experiment than from FM fluctuations. 
Whereas the FM fluctuations give a reasonable Tc within 
a wide range of volumes, and can develop as soon as the 
hep structure is stabilized by pressure, the AFM fluctu- 
ations lead to superconductivity only near the magnetic 
instability. Thus, if AFM fluctuations are responsible for 
the observed Tc it is a coincidence that it appears just 
after the bec-hep transition. 

Two possibilities to explain experiment in terms of 
spin fluctuations can be discussed. First, as supercon- 
ductivity due to FM fluctuations is sensitive to impurity 
scatteringE], it is possible that lattice defects induced 
in the pressure experiment near the structural transition 
can suppress Tc- A decrease of 8-10 K would bring Tc 



into agreement with the measured one within a correct 
pressure range as soon as the hep phase becomes sta- 
ble. Lattice imperfections like dislocations, interstitial 
sites, vacancies etc. may be generated especially near the 
critical volume for the structural structural transitions, 
leading to a gradual onset of Tc- A second possibility 
is that the concentration of imperfections will suppress 
Tc from FM fluctuations completely, leaving the mech- 
anism of AFM-II fluctuations to be responsible for the 
Tc- As was mentioned the Tc and the pressure range 
fit experiment reasonable well, but one must understand 
why AFM fluctuations can resist better to imperfections 
than FM ones. Short wave AFM fluctuations are 'op- 
tical' corresponding to ka = ? in the formula for the 
dispersion. The propagation velocity, ~ dco/dk, is very 
small for such waves, so that they can be confined be- 
tween defects. The scattering with defects is larger for 
dispersive waves with long wave length. The reduction 
of Tc compared to Tc itself is proportional to vp/(£A), 
where vf is the Fermi velocity, I the impurity scatter- 
ing length and A the superconducting energy gap3. Our 
arguments follow if the velocity of the spin excitations 
take the place of Vf- In both of these possibilities it is 
expected that Tc should increase if the concentration of 
defects can be reduced, since superconductivity based on 
FM fluctuations give a larger Tc within a wider range 
of pressures. Further experiments on high impurity sam- 
ples showing superconductivity in a large fraction of the 
volume should be able to distinguish between the two 
possibilities, or if electron-phonon coupling is responsi- 
ble for superconductivity. 
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TABLE I. Paramagnetic and and FM results for hep Fe 
at different lattice constants. DOS at Ef for the param- 
agnetic state, Debye temperature, electron-phonon coupling 
(A ep ), Stoner enhancement, ferro magnetic A s / and exchange 
parameter J. 



a N(E F ) Qd \ep S \ sf J 

a.u. (Ry atom) -1 K mRy / fi% /atom 

4.862 2L3 460 0.37 3.2 0.76 il 

4.766 19.5 550 0.30 2.7 0.54 16 

4.670 17.8 650 0.25 2.4 0.40 20 

4.577 16.2 750 0.21 2.2 0.32 23 

4.485 15.8 860 0.18 2.0 0.24 29 

4.396 13.6 980 0.14 1.8 0.18 44 



TABLE II. Results of spin-polarized calculations for the 
AFM-I configuration at different lattice constants. Sponta- 
neous AFM magnetic moment m, local Stoner enhancement, 
exchange parameter J = d 2 E/dm 2 , and AFM A s /. 



a 
a.u. 


m 

[is 1 atom 


S 


J 

mRy //i B 1 atom 




4.838 


1.09 


9.5 


25 


0.11 


4.814 


<0.02 


30 


2.0 


0.85 


4.766 




13 


3.5 


0.40 


4.670 




8 


8.5 


0.17 


4.577 




4.8 


14 


0.08 


4.485 




3.7 


19 


0.05 


4.396 




3.0 


25 


0.03 



TABLE III. Results of spin-polarized calculations for the 
AFM-II configuration at different lattice constants. Sponta- 
neous AFM magnetic moment m, local Stoner enhancement, 
exchange parameter J = d 2 E/dm 2 , and AFM A s /. 



a m S J A s / 

a.u. /is/atom mRy/z^/atom 

4.766 O90 6 55 O06 

4.670 0.15 10 12 0.18 

4.624 0.09 15 5.5 0.35 

4.600 - 16 4.5 0.34 

4.577 - 14 5 0.27 

4.485 - 8 8.5 0.12 

4.396 - 5 13 0.08 




4.4 4.45 4.5 4.55 4.6 4.65 4.7 4.75 4.8 4.85 



lattice constant (a.u.) 

FIG. 1. Calculated T c from FM (full line), AFM-I (cir- 
cles), and AFM-II (crosses) spin fluctuations in hep Fe as 
function of the lattice constant a. The calculated total en- 
ergy of the bec structure is lower than for hep for a larger than 
about 4.7 a.u.. The difference in pressure between a=4.7 and 
a=4.4 a.u. is about 700 kBar. Superconductivity has been 
observed with Tc below 2 K in a pressure range oij-about 150 
kBar starting near the bec-hep transition pressure^. 
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